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Abstract 



We study the possible regularization of collision solutions for one centre problems with a 
weak singularity. In the case of logarithmic singularities, we consider the method of regular- 
ization via smoothing of the potential. With this technique, we prove that the extended flow, 
where collision solutions are replaced with a transmission trajectory, is continuous, though not 
differentiable, with respect to the initial data. 
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1 Introduction 

In this paper we deal with dynamical systems associated with conservative central forces which 
are singular at the origin. A classical solution does not interact with the singularity of the force, 
i.e., it is a path u € C 2 (T, R 2 \ {0}) which fullfils the initial value problem 



where V(x) G C 2 (R + ,R) is the potential of the force and T denotes the maximal interval of 
existence. As well-known, the two-body problem with an interaction potential V can be reduced 
to a system of this form where u(t) denotes the position of one of the particle with respect to the 
center of mass. Accordingly, we shall term collision the configuration u(t) = 0. As the force field 
diverges at u = 0, collisions are among the main sources of non-completeness of the associated flow. 
This work studies the possible extensions of the flow through the collision that make it continuous 
with respect to the initial conditions. We are concerned with weak singularities of the potential, 
namely logarithms. 

The regularization of total and partial collisions in the N-hody problem is a very classical 
subject and, in the years, different strategies have been developed in order to extend motions 
beyond the singularity [HJ [H [HI [3 EH EH]- Very roughly, these classical methods rely upon 
suitable changes of space-time variables aimed at obtaining a smooth flow, possibly on an extended 
phase space; to this aim, the first step is to determine the asymptotic behaviour of the collision 
solutions and then the phase space is extended either by means of a double covering, or with the 
attachment of a collision manifold. 

In this paper we consider a further, non classical way of extending the flow, related to the tech- 
nique of regularization via smoothing of the potential introduced by Bellettini, Fusco, Gronchi [2]. 



potential. 
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This method consists in smoothing the singular potential and passing to the limit as the smoothing 
parameter e and the angular momentum tend to zero simultaneously but in an independent manner 
(indeed we know that the only collision motions have vanishing angular momentum) . This involves 
an in-depth analysis about the ways the smoothing of the potential coupled with the perturbation 
of initial conditions lead to define a global solution of the singular problem. This technique, when 
successful, has the advantage of being extremely robust with respect to the application of existence 
techniques such as the direct method of the calculus of variation. Let us mention that variational 
methods have been widely exploited in the recent literature in order to obtain selected symmetric 
trajectories for iV-body problems with Kepler potentials ([!]). 

To begin with, we remove the singularity at x — and we denote with V E {x) the smoothed 
function defined as 

V £ (x) = V(Vx 2 +e 2 ) e > 
Then we look at the regularized problem 

r u = w E (|u|) 

P(e) : { (2) 
[ («(0),«(Q)) = (gb,Po) £ K 2 x I 2 . 

Unlike ([T]), the differential equation is no longer singular, so that the initial value problem 
admits a global solution in C°°((— oo, +oo); R 2 ) for every choice of the datum (qo,Po), provided 
VV(a;) is sublinear at infmitjQ. Since we focus on the singularities due to collisions, we fix a ball 
Bq(R) of radius R centered at the origin, where the collision is the only singularity that system 
([1]) can develop and we denote with S(V) C R 2 x R 2 the set of initial conditions (qo,Po) leading to 
collision for the system P with \q \ < R. For every v G S(V) let Up(t) € C 2 (T, R 2 ) be the collision 
solution where T denote the maximal interval of existence such that < R- Denoting with 

u E ^{t) the solution of ^ with initial data v, we investigate the existence of the asymptotic limit of 
the paths u £ ^(t) as (e, v) — > (0, v), its relationship with the collision solution u„(t) of the singular 
system P and the continuity of the limit trajectory with respect to initial data. The definition of 
regularization considered in [2] is the following. 

Definition 1.1. Let V(x) be a singular potential. We say that the problem (Qp is weakly regulariz- 
able via smoothing of the potential in Bq(R) if for every v S S(V) there exist two sequences (ek)k, 
(vk)k tending to and v respectively, such that there exists 

Urn u Ek , Vk = u 

k— >oc 

and the flow 

u v {t) 

is continuous with respect to v . 
In addition we say that 

Definition 1.2. The singular one centre problem ([T]) is strongly regularizable via smoothing of the 
potential if there exists R such that for every v £ S{V) there exists 

lim u ev = uq (3) 

( E) i0-(o,p) ' 

and the flow 

u v {t) 

is continuous with respect to v. 

1 Without any additional assumption on the behaviour of the potential V(x) far away from the origin, a solution 
of system Jy might have singularities other than collisions: for instance solutions could blow up in finite time. 



f u„(t) v£S{V) 
\ uo(t) v g S(V) 



f u„(f) v£S{V) 
\ u (t) v e S(V) 
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In both the definitions we mean that the limit of the regularizing paths u ejll {t) and the continuity 
of the extended flow are held in the ball B (R). 

In the authors prove that in the case of homogeneous potential of degree a, V(x) = 
a > 0, the one-centre problem is weakly rcgularizable via smoothing of the potential if and only if 
a is in the form 

where n is a positive integer or a > 2. On the other hand they show that the homogeneous problem 
is never strongly regularizable via smoothing of the potential. Indeed, a necessary condition in 
order to achieve the uniform limit ^ is that the apsidal angle A9i (u) of a solution of the system (JTJ) 
has to converge to ^ as the angular momentum I tends to zero (see the definition of apsidal angle 
in the next section). This condition is never satisfied by a-homogeneous potentials a > 0, since 
A9i (u) — > 2^ > § as I — > [12] . This explains why the one centre problem with homogeneous 
potential can not be regularized according to definition 11.21 Conversely, when the logarithmic 
potential is considered, it can be proved f[12j) that the limiting apsidal angle do indeed converge 
to 7r/2. Then there is no obstruction and we could expect that the limit §S§ is attained. This 
fact suggests to extend the motion after a collision by reflecting it about the origin. We will show 
that, in this way, not only for the logarithmic potential, but for a larger class V* of potentials, the 
problem is regularizable according to definition 11.21 

The sets of potential functions we will consider in this paper are the following. 

Definition 1.3 (The function set V ). We define V the set of functions V(x) € C°°(]R + ,]R) with 
the properties: 



i. lim V(x) — +oo 
there exists S > such that for every x € (0, S) 

a. v'(x) < o , v"(x) > o 

V'(x) 

iii. the function ^— is decreasing with respect to x 

V"[x) 

and 

iv. V4(o) < -- 

dxV"(x) y ' 2 




The properties iii,iv guarantee the existence of a T > such that 

V'(x) x 

— i-T < n 0, T 

V"(x) ~ 2 V ' ' 



(5) 



Let 



R := min{T, S}. 



(6) 



Definition 1.4 (The set V*). Denote with V* the set of functions V(x) G V with the further 
property 

V(Xx) 

v. lim y7y> = 1 f or ever "!J x > 1 uniformly in every compact K — [1, M\. 
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The set V includes potentials having homogeneous singularities and weaker. For instance the 
logarithmic potential, V(x) = — \og(x), as well as the homogeneous potentials, V(x) — \x\~ a , 
provided a £ (0, 1), belong to V. On the other hand condition v. can be considered as a loga- 
rithmic type property or a zero-homogeneity property: indeed it is never satisfied by homogeneous 
potentials, while the logarithmic potential is a prototype of all the functions satisfying condition 
v.. 

Our main goal is the following: 

Theorem 1. For every V(x) £ V* the one centre problem is regularizable according to definition 
\1.2\ where R is given in © . 

In the particular case of logarithmic potential, V{x) — — log(a;), one has R — +00, therefore 

Corollary 1.1. The logarithmic one central problem is globally regularizable via smoothing of the 
potential according to definition ] 1 . 2[ 

The paper is organized as follows. In section [2] we follow the classical method for dealing with 
central problem based on first integrals and we derive the set S(V) of initial conditions leading 
to the collisions for the unperturbed system. Next, in section [3j given any collision solution u(t), 
we set the initial data v £ S(V) and we define the family of paths u e _ u (t). Section 2] contains the 
proof of the main theorem and the analysis of the regularity of the extended flow. The main part 
of the proof consists in proving the existence of the limit of the path u £tlJ (t) as (e, v) — > (0,£/), 
especially for what that concerns the angular part, (Theorem [5]) . This is the most delicate step, 
for the it involves the uniformity of the limit as (e, v) — > (0, v) 1 and it allows to conclude the strong 
regularizability of the problem. 

It results that the natural extension of the collision solution is the transmission solution, see 
definition 14. 1[ obtained by reflecting the motion through the collision. The regularity of the 
extended flow is carried on in section 14.21 in theorems G2 and H the continuity of the Poincare map 
and Poincare section with respect initial data is achieved. 

Furthermore, in order to have a complete picture of the problem, in section [5] we join a varia- 
tional approach and we analyse the variational properties of the collision paths. 



2 Preliminaries 

For any choice of the potential function V(x) £ C 2 (K + ,M) the one centre problem |T]) is a hamil- 
tonian system and admits two integrals of motion: the energy E and the angular momentum 
I: 

E =^\ii\ 2 -V(\u\) l = uAu 

Since the conservation of the angular momentum implies the motion is planar, in the following I is 
used to denote the modulo of the angular momentum, rather than the vector. The radial symmetry 
of the equation of motion fl} suggests to introduce the polar coordinates in the plane (r, 9) . In 
this setting the quantities E and I are expressed in the form 

E = \i a + ~-V{r) 1 = tH (7) 

We define 

f(r) = 2r 2 (E + V(r)) (8) 

then the relation ([7]) reads as I 2 = f(r) — (rf) 2 and shows that a solution of system ([1]) of energy 
E and angular momentum I exists only for those values of the radial coordinate r > such that 
f(r) >l 2 - 

For fixed values of E and /, we denote with if it exists, the minimum positive value of 
r such that /(?•) = I 2 and f'(r) < 0, R + = +00 otherwise, and we denote with if it exists, 
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(a) Apsidal values 



(b) Apsidal angle 



Figure 1: 



the minimum positive value of r < R+ such that f(r) — I 2 and f'(r) > 0, R- = otherwise, 



see figure 1(a) By definition, it follows that i?_ = for collision solutions and i?+ = +oo 
for unbounded orbits. In any case we term R + and i?_ the maximal and minimal value of the 
angular coordinate and we refer to them as the apsidal values of the orbit. Moreover, following 
the terminology adopted in celestial mechanics, we sometimes refer to i?+ and i?_ respectively as 
the apocentre and the pericentre of the orbit. As it is well known, in case of non collision and 
bounded trajectories, the radial coordinate r(t) oscillates periodically between its extremal values 
R + and i?_ while the angular coordinate 0(t) covers an angle equal to 

R+ i 



A0j(u) = / dr 
R - r^(E + V(r))-l 

between each singular oscillation of r(t). We term the angle A9i(u) the apsidal angle, see figure 



1(b) 



The knowledge of the apsidal values i?+ and i?_ and the value of the apsidal angle is sufficient to 
determine the behaviour of the solution since the whole trajectory is obtained repeating periodically 
the part of path between a point where r(t) is maximum and the following point where r(t) is 
minimal. 

The definition of the apsidal angle extends in a natural way for unbounded and collision so- 
lutions: in the first case the orbits is composed by a single oscillation of the radial coordinate 
from infinity to i?_ and back to infinity and the apsidal angle represents the angle covered by 
the particle coming from infinity to the pericentre and it is obtained replacing R + = +oo in the 
previous relation. On the other hand, if a collision occurs, the apsidal angle denotes the increment 
of the angular coordinate between the apocentre and the collision point. 

In order to characterise the set S(V) of initial data leading to a collision we give the following 
definition. 

Definition 2.1. We say that a potential function V(x), singular in the origin, is of weak type if 

lim x 2 V(x) = 

Otherwise we say that V{x) is a strong type potential. 

A similar classification of singular potentials can be found in a work of Gordon [5] where a 
potential V(x) is said to satisfy a strong force condition at a point xq if V(x) tends to infinity as 
x tends to xq and also there exists a function U(x) with infinitely deep wells at xq, such that 

V{x) > \VU(x)\ 2 

in a neighbourhood of xo- We notice that, among the homogeneous potentials, the set of potentials 
with the property to be of strong type and the ones satisfying the Gordon's strong force condition 
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coincide. Clearly, an a-homogeneous potential is of weak type if and only if a £ (0, 2) and in these 
cases a collision occurs only in zero angular momentum orbits [8] , while if a > 2 a collision solution 
exists also for non-zero values of the angular momentum 2,8. The next proposition extends this 
result. 

Proposition 2.1. IfV(x) €C 2 (R + ,R) is a potential of weak type, a solution u(t) of the dynamical 
system u = VV(|u|) ends into a collision if and only if the angular momentum is zero. 

Proof 

Denote with E and / the energy and the angular momentum of the solution u(t) and let f(r) as 
in ((HI). As mentioned before, a solution exists only for the values of radial coordinate r satisfying 
I 2 < f( r )- Suppose / = 0: since V(x) tends to infinity as x goes to zero, for every value of E there 
exists a neighbourhood of the origin where E + V(r) > then the solution presents a collision. 

Conversely, since V(x) is a weak type potential, it follows that f(r) — > as x — > + thus for 
every value of / ^ there exist a neighbourhood of the origin where I 2 > f(r). Hence the collision 
can not be attained on solutions with non zero angular momentum. □ 

Proposition 2.2. Every V{x) € V is a weak type potential. 

Proof 

From relation (O, by integration, it follows the estimate 

-v'(0 < , Cx > 

for every £ £ (0,T). Therefore, again by integration, we infer 

V{x) <— + C 2 (9) 
x 

and we conclude 

lim x 2 V(x) = . 

□ 

Propositions 12.11 12.21 show that, for any choice of V(x) € V, a solution of the system (fl]) ends 
into a collision if and only if the angular momentum is zero. Therefore the set S(V) of initial 
conditions v = (qo,Po) that make the solutions to be singular consists in those (q~o,po) satisfying 
I = \% A Pol = o. 

3 Setting 

For every fixed V(x) £ V let R be the quantity defined in © and let B$(R) be used to denote 
the ball of radius R around the origin. The properties of the potential class V assure that the 
collision is the only source of singularity for the dynamical system inside Bq(R). As discussed in 
the introduction, given a collision solution u(t) for the one central problem (fl]), our intent is to 
define an extension in Bq(R) for the solution u(t) beyond the collision. 

To this we first have to set the initial conditions v = (qo,Po) for the singular path u(t): we 
denote with P the first positive solution of equation f(r) = 0, P := +oo if such value does not 
exists. We remind that the angular momentum I is zero for collision solution, hence the orbits 
drawn by u(t) is a straight line joining some point in the plane with the origin. 

An alternative occurs: 
Case 1 P < R. 

The collision solution u(t) is bounded in a ball centered in the origin of radius P < R and P 
represents the maximal value of the radial coordinate ( Figure [2]) . Without loss of generality, we 
can set the initial condition v of the collision solution u as 

P = (g,Q), \q\ = P (10) 
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Figure 2: Case 1 



Case2 P>R 

In this case the collision solution is not bounded in Bo(R) and it could also be unbounded. We 
focus our analysis only on the portion of path bounded by R hence we select as initial condition 
for u(t) the couple 

v=(q,jf), \q\=R, \p\ 2 = 2(E + V(R)) (11) 
where the initial velocity p is directed toward the center of attraction ( Figure [3|) . 




Figure 3: Case 2 

In both the cases, for e > let u EtV {t) be a solution of system ^ leading from an initial data 
v = (qo,po). We refer to u StU (t) as regularizing paths in order to underline the purpose to define the 
extension for the singular solution u(t) as a limit of u e ^(t) as (e, v) — > (0, v). The smoothing of the 
potential does not affect the hamiltonian structure of the system, therefore the angular momentum 
I and the energy E € , v = J; |u £j „| 2 — V^(|u e ,„|) are conserved along the solution u e<v (t). Decomposing 
the initial velocity po in terms of the parallel and orthogonal component with respect to qq, 

Pa = Vi + v 2 vi || q , v 2 J- qo (12) 

one has | -^2 1 2 = ]^j2 and E e>v — i | | 2 + |j^j2 — ^(kol) then the condition v — > v is equivalent to 

q -» q, l->0, vx->p 

that, coupled with the condition e — ► 0, implies E £jV — > E. It turns out that also the apsidal 
values of the regularizing paths have to converge to the corresponding ones of u(t), indeed as 
{e,v) — > (0, 9), the pericentre i?_ of the solution u e<v (t) tends to zero while the apocentre i?+ is 
bounded by R and tends to P in case 1, while R+ > R and possibly R+ = +oo in case 2. 

In the following sections we will deal with the existence and the property of the limit for the 
paths u £ _ v {t) as (e,v) — > (0, v). As discussed in [2], the behaviour of the angular coordinate of 
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the regularizing paths plays a fundamental role for the existence of the limit of u Ej „(i) uniformly 
as {e,v) — ► (0, v) rather then for subsequences {ek,Vk) ~* (0, v). For this reason and since we 
focus our analysis only inside the ball Bo(R), we extend the definition of apsidal angle for the 
solution u £ . v (t) as follows: we denote with A9(u £ . u ) the apsidal angle of the path u £ , v {t) in case 
the trajectory is bounded by R, easel 

r R+ i 

A0(u E ,„) - / -^-dr (13) 
Jr. r l r 

otherwise, in case 2, we denote with A6(u e<u ) the angle covered by the path u EM (i) between 
the point where u £ ^ v (t) enter in the ball Bq(R) and the point of minimal distance from the origin 

r R I 

A9(u e , v ) = / -^dr . (14) 
J R _ r 2 r 



4 Proof of Main Theorem and property of the extended flow 

The proof of theorem [T] is composed by two parts: first, in section |4~T1 we prove the existence of 
the limit of the trajectories u ejl/ (t) as (e, v) — > (0, v) and we define the extension of the singular 
solution, then in section [4721 we study the regularity of the extended flow. 

A necessary condition for the existence of the limit is the existence of the limit of the angular 
part of the regularized solutions. Theorem [5] concerns the asymptotic of A6(u e<u ) as (e, v) — ► (0, v): 
to this aim we first prove in lemma 14.41 the L 1 boundness of the integrand in (fT5| and (fTl")) then 
we apply the dominated convergence theorem and pass to the limit under the integral sign. The 
boundness of the integrand is a consequence of a technical estimate stated in the proposition 14.11 
and it is attained for every potential V(x) G V, while the existence of the limit is a consequence 
of proposition 14.31 based on assumption v. . The result we obtain suggest to define the extension 
i*o (i) of the collision solution u(t) beyond the singularity as a transmission solution. 

In order to gain the regularity of the extension, we analyze, in theorems [3] and 01 the continuity 
of the Poincare map and Poincare section of the extended flow in the phase space. 



4.1 The existence of the limit of u e i(t) 



Proposition 4.1. Let V(x) € V and R as in 
that Ve < e and for every < y < f it holds 



Then for every r < R there exists an e such 



p < \— ( r + x ) 

^(ytf) \ £ J x ) ■— _ ^ 



V £ (x)-V £ (f)\ {f 2 ~y 2 ) 



V £ (y)-V £ (f)J (f 2 -^) 



> f 



(15) 



for every y < x < r . 



The proof of the proposition ^. H is split in two parts: we first show that there exists e > such 
that for every e < e the function Fi yP \{e,x) > F^ y ^(0,x) for every 0<y<x<f<S, then we 
show that Ffyf\ (0, x) > f for every 0<y<x<f<R. 

Lemma 4.2. Let V(x) be a function satisfying the properties i.-iii.. Then for every choice of 
< y < f < S there exists e > such that Ve < e and Vie <E (y, f) it holds 



V £ (x)-V £ (f) V(x) - V(f) 



(16) 



Proof 

By means of straightforward calculations and reminding the definition of smoothed potential 
V s (-) = V(y/{-) 2 +e 2 ), the relation (HI is equivalent to 
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VWv 2 + £ 2 ) -V(Vx 2 +£ 2 ) V(y) - V(x) 
V(Vx 2 +e 2 ) - VWf 2 +e 2 ) ~ V[x) - V(f) 

For every s > we define the function J7(s) = V(\/s). Obviously the function U(s) inherit property 
i. and property ii. for every s £ (0, S 2 ), while the relation 

U"( S ) = V"^ 3 ) J_ 
V U'{s) V'Ws) 2Vi 

and property iii. implies that the function ^/s is increasing for every s G (0, S* 2 ). The proof 

U (s) 

of the lemma follows proving that for every choice of < y < f < S 2 there exists e > such 
that Ve < e and Va; £ (y, f) it holds 

U(y + s)~U(x + e) U(y) - U(x) 

U{x + e) - U(r + e) ~ U(x) - U{f) ' 1 ' 

Let < y < f < S 2 be fixed and let the function g(e,x) be defined as g(e,x) := ^jfefOjfepf] i 
x S (y, f). Since g(0) = it's enough to show that 

-^(0,a:)<0 V0< 9 <i<f<5 2 . (18) 
as 

The sign of the derivative is given by the sign of the function 

G (y , P) (x) := U'(y)(u(x) - C/(f)) + U(y)(u'(f) - U'(xj) + U'(x)U(f) ~ U'(r)U(x) 
We observe that G( y ^(y) = G( y ^<)(f) — and 



dx 



(x) = U'(x) (U'(y) - U'(fj) + U"(x) (U{f) - U(y)) (19) 



Since G( y ^(x) is continuous, there exists at least one point x G (y, r) where dG j^' r) (x) — 0; the 
proof of the lemma follows once we prove the inequalities 



' -(x) < for i€ (y,x) 

d< 



4 X ' (20) 



dx 



(x) > for x £ (x,r) 



Indeed, if (|2"U|) hold, the function G( Vir) (x) and, consequently the derivative in (fT5|) are negative 
for every x £ (y,f). To obtain the relations (|20|) we multiply both the sides in (fT9j) for -^a; and 
divide them for U'(x)(U(f) — U(y)). We obtain 

ViG^(x) E/^) (C/'(y) - t/'(f)) 

1 J ' U'(x)(U(f) - U(y)) V f/'(z) + V ((7(f) - C/(y)) 

U"(s) 

By definition of a: it follows N(x) — 0. Moreover, since ^/s . is increasing for every s £ (0, S 2 ) 

and the factor ^\ ~}} is positive, we infer that N(x) is increasing in x. Therefore N(x) < 

{U{r)-U[y)) 

for a; < x and 7V(a;) > otherwise and, since U'(x)(U(r) — U(y)) > 0, inequalities ([2"U|) hold. □ 
Proof of Proposition \4-l\ 
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We fix < y < r < R. For lemma fOl there exists e > such that F( Vi f )) (e, x) > F( y ^(Q,x) for 
every e £ [0,e] and for every x £ [y, f]. Hence it's sufficient to prove 



(r + ai) 



/V(a;)-V(r)\ (f 2 - y 2 ) 



V^(y) - 7(f) / (f 2 -x 2 ) 



- 1 



> f 



Suppose for a moment that relation 



V(x) - V(f) 
V(y) - V(f) 



> 



x\ y 



(21) 



holds for every x £ {y,r). Replacing into Ft f y \{x) we obtain 



(r,y) 

F if , y) (x)>p J '' 



(r + y) 

yj (f + x) 



(f-D 

In order to prove relation (|2ip we rewrite it as 

(V(x)-V(f))_^ (V(y)-V(f)) 



> f 



-x > 



y \/x £ (y,r) 



(r-x) (r - y) 

For x = y the inequality holds; moreover, denoting with N(x) the numerator of the derivative 



(22) 



dx 



{V{x)-V(f)) : 
(f — a;) 



x(f - x)V'{x) + f(V(x) - V(r)) 
(f — a;) 2 



(23) 



one has TV (f) = and ^(a:) = (r-x)(2V'(x)+xV"(x)).For ©, for every a; £ (0, R), £g(x) < 0, 

□ 



hence the derivative in (f2"3")) is positive and relation (f2"2")l holds for every x G (y,r). 
Proposition 4.3. Let V(x) £ V* , then for every p > 1 

Ve(pS) , 



lim 

(<s, s )-(o,ci) K(<5) 



6,s > 



Proof 

We rewrite the above limit in the form 



lim 

(5,6)^(0,0) 



VWS 2 + e 2 ) 

Since p > 1, for every choice of positive values of s and <5, it holds 

p 2 6 2 



1 < 



,2^2 _l_ e 2 

J 2 + e 2 



then, from definition (|1.4[) of V*, replacing A and M with £ 2 + e 2 and p respectively, we infer the 
statement of the proposition. □ 
Let V(x) £ V* and u(t) any collision solution of the system ([T]) with energy E and leading from 
the initial condition D — (cjo,po) in the form (|10p or pip . According with the previous setting, 
section [31 for every sufficiently small e > 0, let U £)! ,(t) be the solution of the regularized system 
^ with initial condition v. Reminding the definition of A8(u EM ) given in (fT3)) . (fTlf we state the 
following theorem. 



Theorem 2. There exists 



and such limit is \. 



lim A9(u e v ) 

(e,iO-KO,P) 
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Proof 

Reminding the definition of R and R + , we define 

(3 = min{i?, R+} 

therefore, for every u Ei i(t), regardless they are bounded or not by R, we write 

" I 



r I 

A9(u e . v ) = / — dr 
J R _ r 2 r 



In order to deal with the convergence of A9(u e ^) we first rewrite the integrand in a different way. 
From the conservation of energy it follows 

r 2 = 2(E £ ^ + V £ (r))-^ 

thus, replacing the radial velocity and extracting the roots R- and j3 from the denominator, we 
infer 



A6{u £yV ) = / dr = \ 

Jr - rJ^iE^ + VJr))-! J* 



1 ff 1 (r -#_)(/? -r) 



(E e>v + 14(r)) - 1 J*- W(r-R.)(P-r)^ ^{E e>v + V e (r)) 1 



By means of the change of variables p = we rewrite the integral in the form 



1 (0-pR-)(p- 1) 
A6(u £ v ) = / — OD2 2 K , — — '-r dp 

K py/W - pR^){p -l)^J^L (E e , + Ve {pR-)) - 1 

Kdp (24) 



1 py/(R- P R.)(p-l) 

We proceed with the proof of the theorem as follows: first we exhibit, for every V(x) <G V, an 
uniform bound for the function K provided e small enough is taken, then we apply the Lebesgue's 
theorem in order to obtain the limit of A9{u £yV ) as (e, v) — > (0, v) and we'll prove that such limit 
exists if V{x) e V*. 

Lemma 4.4. Lef V^rr) e V, i/ien /or e small enough the function K is bounded by a constant in 
its domain. 



Proof 

(/3-pfl_)(p-l) 



Replacing in K the relation E £ , u = \v 2 + \\p — V e {[3) we obtain 



^ %^(hv 2 + h^-VM + Ve(pR-))-l ^ 



We observe that v is the zero in case (3 = R 



Subtracting the energy formula = i|u £i ;| 2 — Kr(w e j) evaluated in i?_ from the same 



evaluated in /3 we infer 

2R 2 _ 1 f3 2 -R 2 _ 



I 2 {V £ {R.)-V e {[3) + \v 2 ) d 2 

that, replaced into if, implies 



K = (26) 

V £ ( / 9fl_)-V E (/3)+^ \ /3 2 -« 2 _ 



v p y V v E («_)-y E (/3)+^- / 
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We note that for every < a < b, the function f(z) = is increasing for positive z. 
The condition ii. in definition (|1.3|) implies that, for every small enough e, the function V e (x) is 
decreasing with respect to x for every x < (3. This yields the relations V e {(3) < V e {pRJ) < V S (R-), 

2 

thus replacing a = V e (pR-) — V e (/3), b — V e (R-) — V s (0) and z = in f(z), it follows 

V e (pR-) -VM + j ^ Ve(pR-) - Ve{P) w A l 



Therefore 

1 



K< 



>R-) ( 2 ( 

p-1 \P \ V s {R-)-V B (p) ) 2 -R± p 2 

and, by means of the substitutions x — pR- and y = R- 

f" 1 \\y) \V e (y)-V e (f3) J f3 2 -x 2 

Finally, applying proposition (|4.1[) . we obtain the uniform bound 

K<[3 \/ p€ {l,JL 



Continue the proof of theorem [H 
The boundness of the functions K \ and K% and the formula 



□ 



i: 



i xJ(x-l)(l-f) 



-.dx = 7T £ > 1 (27) 



implies that all the integral (j2l|) is bounded by n/3. 

In order to apply the Lebesgue dominated convergence theorem we need an uniform bound of 
all the integrands, independently on the values of {e,v) in the neighbourhood of (0, v). We note 
that the singular point moves as e and v change, thus we rewrite A9(u £ iU ) as follows 

£ ' h pV(0-pR-)(p-i) 

1 r Kdp+ [ R ~ — L= =VKdp = h+ 1 2 



P ^{[3- P R-){p~l) JJ3E Py /(/3- P R.)(p-l) 



and 



lim A9(u e v ) = lim I± + lim I 2 

(e,«/)-»(Q,P) ' (e,!/)-»(0,!7) (e,i/)->(0,») 



First we calculate the integral I 2 and we check that it is infinitesimal as (e,is) — > (0, f): for the 
boundness of if, K < j3, 

1 [RT i fsr 1 

Z2 = ~7g / , „ ; = ^=VKdp < I ; = ^=dp 



£ py/(l-p^-)( P -l) J ^= p^(l-p^)(p-l) 
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hence, for ijTTjl . 



h < 7T - 



-dp = tv — 2 arctan 



p J{l-p^)(p-l) 



P-1 



1 - 



pR- 



arctan , 



VP 



/R- 



2 arctan • 



^ 2? R - 



where, in the last passage, the relation arctan(a;) + arctan(i) = -| is used. Passing to the limit, 
reminding that i?_ — > as (e, v) — > (0, v), we infer 



lim I 2 = . 

It remains to prove the existence of the limit for I±. We define F e ^(p) 

1 



the function 



pViP- 1 )^- p r -) i'V^-i 

thus the limit of I\ is equivalent to 



lim 

(e,i/)-K0,P) 



F e , v {p)dp 



We note that every function F etV (p) becomes unbounded only for p approaching p = 1. Moreover, 
for the boundness of -R', it follows that Vp G (1, +oo) 



C 



c 



< 



C" 



then all the functions F e ^ v {p) are dominated by a function F e L 1 ([l,co]). We apply the Lebesgue 
theorem computing the pointwise limit of F eM (p) as {e,v) — > (0,9). 
Using K in the form (|26p we write 



lim 



(e^)^(o,p) py/(p-l){j3- pR_) 



lim — 

(i,e)->(0,0) p 



\ 



ft 2 



and reminding that the convergence (e, v) — » (0, f) implies / — > and, as a consequence, i?_ — > 
we gain 

1 / 



i 



F \ P V E (R-) 1 

Therefore for every V(x) G V*, thank to proposition (|4.3p . we infer 

lim F e i (p) = ; 

(l, e )-(o,o) ' VP/ pVp 5 ^ 

uniformly in e and /. Thus, for the Lebesgue theorem, there exists the limit of A#(u Ejl/ ) and it 
values 

f°° 7T 

lim / F F i(p)dp=—. 
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□ 

Since for (e, v) — ► (0, P) the apsidal angle A9(u E:U ) tends to |, the pointwise limit of the 
sequence of trajectories u e j(t) in the ball Bq(R) is a straight line trajectory that crosses the origin. 

This suggests to extend the collision solution u(t) beyond the singularity replacing symmetri- 
cally the solution itself forward the collision point in the same direction. 

Definition 4.1. Let u(t), t S [0, To), be a collision path, Tq the collision instant. Define the 
transmission solution Uo(t), tE [0,2Tb] as 

f u (t)=u(t) te [0,T ] 

\ u (t)=-u(2T -t) te [T ,2T ] 

In order to complete the proof of theorem Q] it remains to show that, for every t £ [0, 2Tb], 
the sequence {u £jV (t)} pointwise converges to uo(t) as (e,u) — ► (0, v) and that the flow obtained 
replacing the collision solution u with the transmission solution uq is continuous with respect initial 
data. 



4.2 The regularity of the extended flow 

We denote with T = {y = (x, x) £ I 2 x I 2 } the phase space of planar motion and we consider the 
initial value problems defined on T equivalent to systems ((TJ and @ 

{y — f(y) { y = fe(y) 

x(0) = xo e K 2 \{0} , P(e) = I x(Q) = i eM K 
i(0) = ±o e { x(0) = ±o e R 2 

where f(x,x) = (x,VV(\x\)) and f e (x,x) = (x,VV E (\x\ j). 

For every initial data y = [q,p), \q\ < R, leading to collision for the system P(Q), let 
y(t) = (x(t),x(t)) : [0, To) — * T be the corresponding singular solution where T denote the 
collision time and |x(£)| < R for every t € [0, To). We extend y(t) according to the definition 14.11 
defining y (t) = (x (t), x (t)) as 

yo(t) = y(t) ie[0,T o ) 

yo{t >-\ x (t) = x(2T -t) * e ( T o,2T ) 

The extension for the collision solutions allows to define the flow beyond the singularity, therefore 
for e > we denote with t) : T x K + — > JF the flow associated to the system P(e). 

Our aim is to study the continuity of <f> e (y,t) with respect to initial data y € T and e as 
(y,e) — * (y, 0); to this end we consider the Poincare map $r(j/, e) defined as the solution at time 
T of the system P(e) with initial value ?/ and we show that 

lim $ r (j/,e) = $t(j/,0) 

for every T ^= Tq. 

We note that the continuity of the Poincare map implies the proof of theorem [T] 

Remark 4.1. We can not expect the continuity of the Poincare map in Tq because, even if the 
configurations XT (y,s) would converge to XT Q {y,0), the limit can not be attained by the sequence 
±T (y,s), since XT Q (y,0) is unbounded. 

If T < To no problem arises, indeed the above limit comes from the classical theorem of 
continuity with respect initial data of ordinary differential equations. 

On the other hand, for T > Tq the continuity of the Poincare map is stated in the following 
theorem. 
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Theorem 3. Let V(x) G V* and suppose y = (q,p), \q\ < R be an initial condition leading to 
collision for the system P(0) at time Tq. Then 

lim $ T (y,e) = *t(j/,0) 

for every T > Tq such that x(T) =/= . 

In the proof of the theorem we will need the following classical lemma 

Lemma 4.5. Let H(r,p) and Hq(p) be real continuous functions with respect to a set of parameters 
p, and suppose H be strictly increasing as function of r. Then for every T the function r = r(T,p), 
implicit solution of equation 

T = H (p) + H(r(T,p),p) 

is continuous as function of p. 
Proof of theorem ([3]) 

As usual we set (r, 9) the polar coordinates of the plane then a point y = (x, x) in the phase space 
is replaced by 

x = (r cos 9, r sin 9) 

x = (f cos 9 — r9 sin 9, f sin 9 + r9 cos 9) 

From the definition of ^t{v, s) remain well defined the functions and rx(y, s), 9r(y, s), tt(Ui s), 0t(?/, e) 
denoting, respectively, the values of the radial and angular coordinate and their velocity at time 
T 7^ To of a solution of system P(e) with initial data y. The continuity of $7^(2/, e) is equivalent 
to the continuity of each of the previous functions. 

Let y = (f, 9, f, 9) be an initial data in the phase space leading to collision with nonzero initial 
velocity, f < 0, and denote with E = if 2 — V{f) the energy of the collision solution. 

A point y £ T, y = (ro, 9q, t-q, #o)j tends to y if it holds 

\r - f\ , 1 O -d\^0 mod 2tt 

We start showing the continuity of the function ry(?/,e) as y — > y and e — > 0. The value of 
rr(y, e) is governed by the differential equation r = \J 2(E + V e (r)) — thus rx{y, e) is a function 

of the initial position ro, the couple E, I and the parameter e. Define To(E,l,ro,e) as the time 
necessary to the solution r(y,t) to reach the minimal value i?_ = R-(E,l), then 

T = - :dp + / dp 

Jn- ^2{E + V £ {p))~^ Jr - ^2{E + V £ {p))-$ 



= T {E,l,r ,e)+T(r T (y,e),E,l) 



Claim The function To and T are continuous respect to ro, E, I, s. 

Suppose for the moment that the claim is true, since the function T is strictly increasing with 
respect to fr(y, e), for lemma (|4.5|) . the function rr(2/,e) is continuous with respect to the set of 
parameters E, I, ro, e. Therefore 

lim r T (y,e) = r T (y) . 
£—0 

The continuity of 9t(v,£) is equivalent to the continuity of A9r(y,£) '■= 0t(u,£) — &o- F° r the 
definition of transmission solution A9t(v) = 7T, while, in theorem we showed that the apsidal 
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angle associated to a solution of the perturbed system tends to -| whenever the initial data tends 
to a colliding one and the parameter e tends to zero. Thus we gain 

lim 9 T (y,e) = 6 T (y) 

The continuity of rT(y,e) and 8t{u,£) follows immediately by the continuity of rT(y,e) and 
relations 



r T (y, e) = 2 X E + V £ (r T (y, !* x 2 , T (y,e) 



2r T (y,e) 27 ' r T (y,e) 2 

□ 

Proof of the claim 

Denoting with p and p the sets of parameters p = (E, I, r , e) and p = (E, 0, F, 0), we have to show 
that 

lim T (p) =T (p) . 

We want to apply the dominated convergence theorem and pass the limit under the integral sign 
in 

r° i 

lim / rfp 



To this aim we first exhibit an L 1 bound for the integrand function. We observe that the integrand 
become unbounded only for p approaching to R- since we have chosen the initial value r < and, 
as a consequence, we can suppose r < 0. 
Let [3 € (R-,ro) and write 



dp = i ^ dp + / 1 rfp 



H- ^ 2 (E + y e (p))-£ •/*- /2(£ + K(p))-^ ^ j 2 {E + V s {p))-% 

= 1 + 11 

Since for every p e [/?, r ] the integrand is bounded, it follows II < Ci(r — [3). In order to show 
the L 1 boundness of the first integrand, using the definition of the energy integral, we rewrite I in 
the form 

rP i rP i 

1=1 dp = I — dp 

jR - h - Jrt+l>(^f)+2(V £ (p)-VAro)) 



Again from the definition of energy, evaluated in ro and in the pericentre i?_ 

1 .o 1 P_ „, , H 2 

2 

it follows 



i 2 
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Figure 4: Poincare section 



and, replacing the last relation in the integrand, we obtain 







«- ./#2 



1 



fli(^o 2 -P 2 ) 



p 2 (^--Ri) 



+ £$Eg)2(V B (ro) ~ ^(^-)) + 2 (^(P) - K(r )) 



rip 



/3 



1 1 



< / = i = — 4° 



V Y ( r 0~' R -) Vi(p)-Vi(r ) 



For proposition 14 . 1 1 there exists a positive constant K such that, for every e small enough, 
1 i? 2 _ (rl-p 2 ) V e (R-)-V e (r ) ^ Kr {p - R-) 



> 



p 2 (r 2 - R 2 _) V e (p) - V £ (r ) r_(r + p) + Kr (p - R-) 
thus 

I < CI —j^^=—j^^^=^^^= < Ci I dp < oo 

Jr- VP - R - j2(V(/3) - V(r )) Jr - yP ~ R - 

In order to obtain an uniform bound for every choice of R-, we perform the variable change 
z = p — R_ and conclude 

I < C 2 -^dz < oo 



We now apply the Lebesgue theorem and we obtain 

lirn / — =dp = / — dp 

v^Jh- ^2(E + V £ (p))~^ Jo V%E + V(pJ) 

thus the function To and, for the same reason, the function T are continuous with respect to the 
parameter p. 

□ 

Let again y be an initial condition leading to collision for the system P(0), To the collision 
instant and, for a choice of T > To, we denote y± = <!>(?/, T) according to the definition I4.ll of 
extension of collision solution. Given £ an hyperplane in the phase space passing through y\ 
and transversally to the flow, we show that for every initial data y near y, there exists a time 
t = r(y) when the trajectory $(y,t) intersects the hyperplane S. As the data y changes, the trace 
S(y) = ^{y, i~(y)) drawn on X is called the Poincare section of the flow on E, see figure HI In the 
next theorem we prove the continuity of the map r(y) and the continuity of the Poincare section 
in a neighbourhood of y. 
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Theorem 4. Let V(x) G V*. Then there exists a 8 > and a continuous function r(y) defined in 
a 8 -neighbourhood of y, Ng(y), such that r(y) = T and 

%r(!/))6S 

Moreover the Poincare section is continuous in N$(y). 

Proof 

Let F a vector in the phase space such that F ■ f(yi) > and consider the hyperplane 

S = {y ; (y - yl) ■ F = 0} 
By definition of y\ it holds (<&(?/, T) — yl) • F = and, since 

t = f( Vl ) -F>0 
there exists an £ > such that 

($(y,T-£)-yl)-F<0 
my,T + Z)-yl)-F>0 

For theorem © and for the sign permanence theorem it follows that there exists a 5 > such that 
for every \y — y\ < 8 

($(y,T-£)-yl)-F<0 

($(y,T + O-yi)-i ;i >0 

For every fixed y the function ($(y, <) — yl) • F is increasing in i: indeed 

^($(y,i)-yl)-F = /(<%,*)) -F 

thus the continuity of the vector space f(y) out of collision points and the continuity of the orbit 
<&(y,i) with respect to both the variables imply that f(Q(y,t)) ■ F > 0. Therefore for every 
y E N$(y) there exists a time r(y) such that $(y,r(y)) € E and t(j/) = T. In order to prove the 
continuity of r(y) we define 

H(y,T) = my,T)- yi )-F 

then r(y) is the implicit solution of H(y, r(y)) = 0. Since H is continuous in y and it is continuous 
and increasing with respect to r, for lemma (|4.5|) the function r(y) is continuous. Moreover, for 
composition of continuous functions we infer the continuity of the Poincare section. 

□ 

5 Variational property of the collision solutions 

In this section we join a variational approach that consists in seeking solutions of the system (JTJ) 
as critical points of the action functional 

A(u) = J C(u,u)dt 

where 

C{u,u)= l -\u(t)f + V{\u(t)\ 



-my lt )-yl)-F 
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is the lagrangian function associated to the equation of motion. This method is well known in 
the literature and it has been extensively exploited in order to find periodic solutions for the TV- 
body problem, see for instance [TJ [H [5] and the references therein. Besides the discussion about 
the existence of minimal paths, an interesting question is whether the collisions are avoided by 
the minimal paths even in presence of weak potentials whose contribution one could expect to be 
negligible by a variational point of view. 

In the following theorem we prove that, despite of the weakness of the singularity of every 
potential V[x) S V*, a minimal path for the action functional can not have a collision in the 
interior of its domain. 

Theorem 5. For every V(x) € V* , let uq : [— T, T] — > R 2 be the extension of a collision solution of 
system (JTJ) according with the definition \4-l\ of transmission solution. Then uo(t) is not a minimal 
path for the action functional A\\, where 

A = {«(*) : u e L 2 ([-T, T]), u{-T) = u Q (-T), u{T) = u Q {T)} 

denotes the set of paths joining the end points of uq . 
Proof 

In order to prove that uo(t) is not the minima for the action A among all the paths in A, we 
perform a variation on the trajectory uq that removes the collision and makes the action decrease. 
Let u\(t) = uo(t) + v s (t) where v 5 (t) is the standard variation 

v{t)= {m) s ^<iti<T 

directed orthogonally to uo(t). 

Let us compute the difference AA = A(uq) — A(ui) and show that for every S sufficiently small 
A.4 is positive: this means that uo(t) is not a local minimum of the action functional. 



(c(u ) - c( Ul ))dt = J T ^ Qki 2 + v(M)) - Ql«il 2 + nM)) dt 

1 (\ii Q \ 2 - \ Ul \ 2 ) dt+ [ V(\u \)-V(\ Ul \)dt 



AA 



Ji Zl J 

= AIC + AV 

We study separately the kinetic and the potential contribute. 
Since the variation v s is directed orthogonally to uq, we gain 

,. Ml 2_ f \Mt)\ 2 1*1 <T! 

|UlWI -\ \Mt)\ 2 + w^kr < |t| < t 



then 



T i 



We show that, for every 5 small enough, the contribution of the potential part is bigger than the 
penalising contribution, due to the kinetic part. 

AV = 2 ( V{\u \) - V{\ui\) dt > 2 / 1 V(\u \) - V(\ui\) dt 



• I 1 V(\u \)-V(^/u 2 + S 2 )dt 
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For every t fixed 

fl d 



hence 



AV>2 / v v ' S 2 d£_dt = S 2 / v ; s ; d£dt 

lo 2^2+^2 J J ^uJ+W 







= S 2 f 1 f s (t)dt = S 2 Rs . (29) 
Jo 

For 5 sufficiently small the functions fs(t) are positive and for Fatou's Lemma 



Wo Jo v^f+C^ "^o |«o(*)| 



>c/ -y'(|« (*)|) <te = Climr(t) = +oo (30) 

The last relation holds since the radial differential equation r(t) = —V'(r(t)) is satisfied, being the 
collision solution uq a radial solution. Combining (|28l) , ([29]) and (|30|) we conclude 

AA = A/C + AV = <5 2 (- + >0 

for every 5 small enough. □ 
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